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At small layer separations, the ground state of a v = 1 
bilayer quantum Hall system exhibits spontaneous interlayer 
phase coherence. The evolution of this state with increasing 
layer separation d has been a matter of controversy. We re- 
port on small system exact diagonalization calculations which 
suggest that a single phase transition, likely of first order, 
separates incompressible states with strong interlayer corre- 
lations from compressible states with weak interlayer correla- 
tions. We find a dependence of the phase boundary on d and 
interlayer tunneling amplitude that is in very good agreement 
with recent experiments. 

The ground state of a two-dimensional monolayer elec- 
tron system at Landau level filling factor v = 1 is a sin- 
gle Slater determinant described exactly by Hartree-Fock 
theory and is a strong ferromagnet with a large gap E g 
for charged excitations 0J^]. This elementary property 
has rich and interesting consequences for the physics of 
bilayer quantum Hall systems at the same total ^, conse- 
quences that are readily appreciated when a pseudospin 
language ||[l) is used to describe the layer degree of free- 
dom. When the layer separation d goes to zero, interac- 
tions between layers are identical to interactions within 
layers. The pseudospin bilayer Hamiltonian is then iden- 
tical to the single layer Hamiltonian with spin and its 
ground state has pseudospin order and a finite charge 
gap. For infinite layer separation, on the other hand, the 
bilayer system reduces to two disordered, compressible, 
uncorrelated v — 1/2 systems. This Letter concerns the 
evolution of bilayer quantum Hall systems between these 
two extremes. 

For small layer separations the difference bet- 
ween interlayer and intralayer interactions breaks the 
pseudospin-invariance of the Hamiltonian, resulting in an 
incompressible easy-plane pseudospin ferromagnet. In 
physical terms the pseudospin order represents sponta- 
neous phase coherence between the electron layers. Sev- 
eral scenarios have been proposed for the evolution of 
the ground state as the layer separation increases fur- 
ther. In Hartree-Fock theory Q, spontaneous interlayer 
coherence is lost if the layer separation exceeds a critical 
value, and the ground state at large separations consists 
of weakly correlated Wigner crystals. While possibly in- 
structive, this picture is known to be incorrect at large 
d since half-filled Landau levels do hot have crystalline 
ground states. Working in the other direction, Bones- 
teel et al. started || from the composite fermion theory 
of isolated compressible v = 1/2 layers, and concluded 
that coupling would lead to pairing between compos- 



ite fcrmions in opposite layers and also, implicitly, to 
a charge gap. Since the pseudospin ferromagnet pos- 
sesses particle-hole rather than particle-particle pairing, 
however, this picture still implies that at least one phase 
transition occurs as a function of layer separation. In a 
numerical diagonalisation study He et al. B predicted, on 
the basis of the system parameter dependence of overlaps 
between exact groundstates and two different variational 
wavefunctions, the existence of two distinct incompress- 
ible states separated by a region of compressible states. 

Experiments, on the other hand, have tended to 
be consistent with the proposal || that a single- 
phase transition from an incompressible to a compress- 
ible states occurs with increasing layer separation at any 
value of the interlayer tunneling amplitdue. Very re- 
cently, in an intriguing new experiment by Spielman et 
al. 0. the tunneling conductance across the layers was 
studied in a sample with extremely small tunneling am- 
plitude. When the ratio of layer separation and magnetic 
length was lowered (at fixed filling factor) below a critical 
value, the conductance showed a very pronounced peak 
around zero bias voltage between the layers, that pro- 
vides direct evidence || for spontaneous interlayer phase 
coherence. This is because in the coherent state, the 
layer index of each electron is uncertain and only in this 
case can tunneling leave the system in or near its ground 
state so that there is no orthogonality catastrophe and 
tunneling can occur at zero voltage. 

Since the critical layer separation found by Spielman 
et al. is close to the one obtained earlier by Murphy 
et al. for the onset of the quantum Hall effect 0, ex- 
periment demonstrates that for vanishing tunneling am- 
plitude the phase transitions at which pseudospin order 
and the charge gap are lost are either closely spaced or 
coincident. 

In this Letter we report on small system exact diag- 
onalization calculations which strongly suggest that bi- 
layer quantum Hall systems have a single phase transi- 
tion, likely of first order, as a function of d. Our crit- 
ical layer separation is in very good quantitative agree- 
ment with the value measured in Ref. In the light 
of the experimental results mentioned above, our calcu- 
lations imply that the charge gap disappears and long- 
range phase coherence simultaneously drops sharply to 
near zero at the phase transition. This result is not en- 
tirely unexpected since a simple Landau- Ginzburg anal- 
ysis indicates that the two order parameters could not 
vanish simultaneously without fine-tuning, if the transi- 
tion were continuous. Also the mean-field theory energy 
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gap is proportional to the pseudospin order parameter, 
suggesting that these two orders are mutually reinforcing 
and that a first order transition is therefore likely. Fi- 
nally, we note that, experimentally, the charge gap phase 
transition is sharp even at finite tunneling between the 
layers. Since tunneling produces a pseudomagnetic field 
which couples to the pseudospin order parameter, this is 
an unusual magnetic transition which does not involve 
symmetry breaking, a fact which lends further weight to 
the suggestion that the transition is first order. 

We analyse bilayer quantum Hall systems numerically 
by means of exact diagonalisations of finite systems un- 
sing the spherical geometry. We have verified numerically 
that the ground state and low-lying excitations are fully 
spin-polarized and neglect the spin degree of freedom in 
the present discussion. The Hamiltonian is given by 

H = Hip + Hcoui , (1) 

where Ticoui represents the usual Coulomb interaction 
within and between layers, and the single-particle Hamil- 
tonian Hip is given by 

Hlp = - jE^^'+Kc'lv,"" ( 2 ) 

TCI 

We concentrate here on the tunneling amplitude (A t ) 
tuned phase transition, although bias voltage (A„) de- 
pendence is also interesting and often experimentally 
more convenient. [i,[t! 6 {+? ~ } run over the layer 
(or pseudospin) indices and a summation convention is 
implicit; f are the pseudospin Pauli matrices, m 6 
{— N$/2, . . . , N<p/2} is the z-projection of the orbital an- 
gular momentum of each electron in the lowest Landau 
level, where is the number of flux quanta penetrating 
the sphere. In the following we denote the pseudospin op- 
erators by T = (1/2) £ m c+ 

m T fj..fj.' c fj.' .m- The interlayer 
separation d is measured in units of the magnetic length 
Ib = yUcfeB, and all energies are given in units of the 
Coulomb energy scale e 2 /e^. We consider the case of 
zero well width to enable comparison with most previ- 
ous theoretical investigations (h]-[]1|), and also systems 
consisting of two rectangular wells of finite width w |l4| ] 
whose ratio to the center-to-center layer separation d is 
w/d = 0.65. This value corresponds to the sample used 
in Ref. §. 

We consider systems with an even electron number iV 
which leads to a nondegenerate spatially homogeneous 
ground state with total angular momentum L = 0. For 
simplicity, let us first examine the case of vanishing bias 
voltage, where both (T v ) and (T z ) are strictly zero. 

Figure |l| shows the interlayer phase coherence as mea- 
sured by the expectation value (T x ) along with the fluc- 
tuation AT X = \J{T x2 ) - (T X )(T X ) as a function of the 
tunneling gap for a system of twelve electrons, a layer 
separation of d = 1.80, and zero well width. At A t = 0, 
(T x ) is necessarily zero in a finite system. With increas- 
ing tunneling gap, (T x ) grows rapidly reaching a inflec- 
tion point with a very steep tangent. The differential 



pseudospin susceptibility, \ = (l/N){T x )/dA t , is plot- 
ted in the inset and shows a very pronounced peak. In 
the immediate vicinity of this peak, the pseudospin fluc- 
tuation AT X has also a pronounced maximum. In figure 
U the x is plotted for different numbers of electrons. 
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FIG. 1. The pseudospin expectation value (T x ) and the 
fluctuation AT X as a function of the tunneling gap A t . The 
derivative d{T x )/dA t (measured in units of 1/ (e 2 / eIb)) is 
shown in the inset. 



The rapid growth with increasing system size of the 
peak in this generically intensive quantity is strong ev- 
idence for a ground state phase transition. Analogous 
findings are obtained for the peak in the pseudospin fluc- 
tuation. Thus, the peaks in the susceptibility of the pseu- 
dospin and its fluctuation grow very rapidly with increas- 
ing system size and signal a quantum phase transition at 
the critical value of the tunneling gap. At large tunneling 
the system pseudospin magnetisation is close to its max- 
imum value, while at small (but also finite) tunneling the 
the system is disordered and the pseudospin magnetisa- 
tion is strongly reduced by interactions. 
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FIG. 2. The pseudospin susceptibility \ f° r different sys- 
tem sizes as a function of the tunneling gap. The rapidly 
growing peak indicates a quantum phase transition. 



The two peaks described above occur at extremely 
nearby values of A 4 at a given layer separation d, and 
we consider the very tiny differences in their location as 
a finite-size effect. To estimate the phase diagram of the 
system we place the phase boundary at the maximum of 
the quantum fluctuations AT X . 



5.0 



4.0 - 



3.0 - 



2.0 - 



0.0 







N=6 


N=6 










,,N=10 
N=8 yVy 






^=12 








w/d=0.00 




w/d=0.65 



0.00 0.02 0.04 0.06 0.08 0.02 0.04 0.06 0.08 
\ \ 

FIG. 3. Phase boundaries for different system sizes N and 
ratios of well width to layer separation. 



Figure || shows the resulting phase boundaries for dif- 
ferent system sizes and both cases of well width. At small 
layer separation the system is in the ordered phase and 
the fluctuation peak occurs exactly at A t = 0. At a crit- 
ical layer separation d c (A t — 0,N) the phase boundary 
moves out rapidly to finite values of A t and intersects 
the axis A t — with an almost horizontal tangent. This 



is in qualitative agreement with earlier experimental Q| 
and theoretical || estimates of the phase diagram. 
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FIG. 4. The critical layer separation d c (A t = 0,iV) (filled 
symbols) at vanishing tunneling as a function of the system 
size N for both cases of well width. The lines are finite-size 
fits to the data with a shift exponent of A = 5.0. 



The critical values d c (At = 0, ./V) form a rapidly con- 
verging data sequence and are plotted in figure ^|. These 
finite-size data are accurately and consistently described 
by an ansatz of the form d c (N) = a + f3N~ x with two fit 
parameters a — d c (N — oo), (3, and a shift exponent A. 
The best fits to both sets of data are obtained for A — 
5.0 ± 0.2 leading to a value of d c (N = oo) = 1.30 ± 0.03 
for zero well width, and d c (N = oo) = 1.81 ± 0.03 for 
w/d = 0.65. The latter value is in excellent agreement 
with the results of Ref. § , where the onset of the tun- 
neling conductance peak is observed at a layer separa- 
tion of d — 1.83. Thus, our numerical results clearly 
indicate that the findings of the above tunneling exper- 
iments are the signature of a quantum phase transition. 
The very large value of A seems inconsistent with a di- 
verging correlation length and suggests the transition is 
first order. A first order phase transition would explain 
the apparent coincidence of the appearance of sponta- 
neous phase coherence and the quantum Hall effect in 
experiment §,§ We note that our result for the critical 
layer separation at vanishing tunneling gap agrees rea- 
sonably, at zero well width, with the point at which the 
uniform density phase coherent state first becomes un- 
stable in the Hartree-Fock approximation ||. At larger 
w, however, the Hartree-Fock estimates clearly deviate 
from the exact diagonalisation result. 

In order to further investigate the order of the quantum 
phase transition, we introduce the ratio 
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where the subscript N refers to the system size. As we 
discuss below, this type of ratio should prove to be a pow- 
erful general tool in the analysis of any quantum phase 
transition. In classical physics this ratio of fluctuation 
to susceptibility is equal to the thermal energy fc^T and 
vanishes at T = 0. The classical relationship does not ap- 
ply here since the Hamiltonian fails to commute with its 
derivative with respect to A t . There is, however, a closely 
related zero-temperature relationship with the typical ex- 
citation energy ojn taking over the role of temperature. 
The fluctuation can be written as 



correlated states also occur with increasing bias poten- 
tial. We predict measureable anomalies in the double- 
layer system capacitance at bias tuned phase transitions. 
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where the sum is performed over all excited states, while 
for the derivative of the pseudospin magnetisation one 
finds from linear response theory 



d(T x ) 
dAt 
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E n — Eq 



(5) 



From these equations we see that ll>n is a harmonic av- 
erage of excitation energies [E n — Eq), weighted by the 
factors |(n|T x |0)| 2 . In particular, ujm has a vanishing 
thermodynamic limit if at least one state with a nonvan- 
ishing matrix element (n|T x |0) has an excitation energy 
(E n — Eq) which extrapolates to zero for N — > oo. Thus, 
equation (^) defines a characteristic energy scale of the 
system at the phase boundary. The operator T x natu- 
rally enters this expression since it couples to a control 
parameter driving the phase transition. 

For a continuous phase transition one would clearly ex- 
pect u)n to vanish at the phase boundary for an infinite 
system, while a finite limit limjv^oo wn is indicative of a 
finite energy scale, i.e. a first order transition. From our 
finite-size data for w/v (evaluated at vanishing tunneling 
and d = d c (N)) we conclude that this quantity extrapo- 
lates for N — > oo to a rather substantial non-zero value 
of order 0.05e 2 /eZ_B ~ 5K for both values of w considered 
here. Along with the arguments and experimental find- 
ings given so far, this result strongly suggests that the 
bilayer quantum Hall system at filling factor v = 1 un- 
dergoes a single first order phase transition as a function 
of the ratio of layer separation and magnetic length at all 
values of the tunneling amplitude. The phase boundary 
separates a phase with strong interlayer correlation (and 
a finite gap for charged excitations) from a phase with 
weak interlayer correlations and vanishing E g . 

Finally we comment briefly on the influence of a bias 
voltage between the layers. When applying a bias voltage 
to the system the vector (T) is tilted out of the icy-plane 
with a finite ^-component. In this case we find numeri- 
cally that the quantum phase transition is again signaled 
by the longitudinal fluctuation of the pseudospin mag- 
netisation and its susceptibility, and all results concern- 
ing the phase boundary are qualitatively the same. First 
order phase transitions from stronly correlated to weakly 
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